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The duality relation between one-loop integrals and phase-space integrals, developed in a 
previous work, is extended to higher-order loops. The duality relation is realized by a modi- 
fication of the customary +zO prescription of the Feynman propagators, which compensates 

^ I for the absence of the multiple-cut contributions that appear in the Feynman tree theorem. 

We rederive the duality theorem at one-loop order in a form that is more suitable for its 
iterative extension to higher-loop orders. We explicitly show its application to two- and 

S^ . three-loop scalar master integrals, and we discuss the structure of the occurring cuts and the 

H [ ensuing results in detail. 
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1 Introduction 



With the start of the LHC, the physics of elementary particles enters a new era, opening a powerful 
window to discover the Higgs boson and to explore new interactions beyond the Standard Model (SM) 
at the TeV energy scale. Precision theoretical predictions for background and signal multi-particle hard 
scattering processes, in the SM and beyond, are mandatory for the phenomenological interpretation of 
experimental data, and thus to achieve a successful exploitation of the LHC physics programme. 

While leading-order (LO) predictions of multi-particle processes at hadron colliders in perturbative 
Quantum Chromodynamics (pQCD) provide, in general, a rather poor description of experimental data, 
next-to-leading order (NLO) is the first order at which normalizations, and in some cases the shapes, of 
cross sections can be considered reliable, yj]. Next-to-next-to leading order (NNLO), besides improv- 
ing the determination of normalizations and shapes, is also generally accepted to provide the first serious 
estimate of the theoretical uncertainty in pQCD. Despite the relatively smaller coupling, electroweak 
(EW) radiative NLO corrections might also be sizable at the LHC, [|2l|3l. 

Computing higher-order corrections in Quantum Field Theories, in particular in QCD or in the EW 
sector of the SM, is highly challenging and substantially demanding as the complexity increases with the 
number of external particles, and the order in perturbation theory at which the hard scattering process 
must be calculated in order to match the experimental precision. In the recent years, important efforts 
have been devoted to developing efficient methods able to boost forward the calculational capability both 
at the multi-leg and the multi-loop frontier. Today, 2—7-4 processes at NLO, either from Unitarity based 
methods, [l4l[5l[6l|, or from a more traditional Feynman diagrammatic approach, [7J, are affordable and 
are even becoming standardized. There has also been a lot of progress concerning NNLO calculations 

In Ref. lfT2l . a duality relation between one-loop integrals and phase-space integrals has been demon- 
strated. The duality relation is suitable to numerically compute, [fT3l . multi-leg one-loop cross sec- 
tions in perturbative field theories (local and unitary). It has analogies with the Feynman tree theorem 
(FTT), lfT4l[T5l , but involves only single cuts of the one-loop Feynman diagrams. The duality theorem 
requires to properly regularize propagators by a complex Lorentz-covariant prescription, which is dif- 
ferent from the customary -\-iO prescription of the Feynman propagators. The main consequence of this 
new prescription is that the multiple cuts appearing in the FTT are avoided. 

The computation of cross sections at NLO (or NNLO) requires the separate evaluation of real and 
virtual radiative corrections. Real (virtual) radiative corrections are given by multi-leg tree-level (loop) 
matrix elements to be integrated over the multi-particle phase-space of the physical process. The loop- 
tree duality at one-loop presented in Ref. lfT2l . as well as other methods relating one-loop and phase- 
space integrals [|T6l[T7llT8l . have the attractive feature that they recast the virtual radiative corrections in 
a form that closely parallels the contribution of the real radiative corrections. This close correspondence 
can help to directly combine real and virtual contributions to NLO cross sections. In this paper, we 
extend the loop-tree duality theorem derived in Ref. [fT2ll to higher-order loops, as a first attempt towards 
extending the duality method to the computation of cross sections at NNLO or even higher orders. 
Preliminary results were presented in Ref. [|19ll . 

The outline of the paper is as follows: In Section [21 we reconsider the tree-loop duality theorem 
at one-loop. This involves the definition of dual propagators in addition to the Feynman, advanced 
and retarded propagators commonly known. Afterwards, we reformulate the duality theorem in a way 
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which is more appropriate for extending it to higher loop orders. This is done by providing functions of 
propagators of complete sets of momenta, corresponding to the internal lines of the diagram. As for their 
single-momenta analogues, it is likewise possible to infer relations amongst these functions. The proof 
of the main relation which is crucial for the extension to higher loop orders, is given in the Appendix lAl 
In Section [3l we then provide a duality theorem for the two-loop master diagram with A^ external legs 
by using the previously defined relations and iteratively applying the duality theorem to the occurring 
loops. We also discuss some subtleties involved, as well as the structure of the result and occurring cuts. 
Furthermore, we derive a two-loop representation of the Feynman tree theorem. Having set up the basic 
method in this way, we will continue with the four basic master topologies at three loops in Section IH 
and show that the method is indeed extendible to even higher loop orders in an iterative manner. We 
end this section with a brief comment on the extension of the duality theorem at the amplitude level. In 
Section [5l we conclude and provide an outlook. In the Appendix |Bl we provide a simple example of a 
two-loop scalar integral calculated from its dual representation. 



2 Duality relation at one loop 

In this section, we provide the basic quantities, definitions and relations used in the rest of the paper, and 
sketch the steps for the derivation of the tree-loop duality theorem at one-loop, as presented in Ref. [[T2ll . 
We will also rederive this tree-loop duality theorem in a form which is more suitable for its extension to 
higher orders by introducing the main formulas used for the iterative application of the duality. 

Let us start by considering a general one-loop A^-leg diagram, as shown in Fig. [U which is repre- 
sented by the scalar integral: 

N 

The four-momenta of the external legs are denoted pi, i G {1,2,. ..A^}. All are taken as outgoing and 
ordered clockwise. The loop momentum is £i, which flows anti-clockwise. The momenta of the internal 
lines Qi, are defined as 

qi = ii+Pi,i, ieai = {l,2,...iV}. (2) 

As is commonly used, we define Pij = Pi + pj+i + . . . + pj. Momentum conservation is equivalent to 
Pi,N = 0. We use dimensionally regularized integrals with the number of space-time dimensions equal 
to d, and introduce the following shorthand notation: 



d'^ 



'I 



(27r)' 



(3) 



The space-time coordinates of any momentum k^ are denoted as k^ = (ko, k), where ko is the energy 
(time) component of fc^j. The Feynman propagators GpiQi) in Eq. ([T]) have real internal masses nii: 

GF{qi) = ^T— - . (4) 

qf -mf + zO 

The derivation of the duality theorem is exactly the same regardless of the internal lines being massive or 
massless (nii = 0), as long as the masses are real. Non-vanishing internal real masses only account for a 




Figure 1 : Momentum configuration of the one-loop N -point scalar integral. 

displacement of the poles of the propagators along the real axis, which does not change the derivation of 
the duality theorem, as will become obvious in the following. Moreover, they do not alter the relationship 
between Feynman, advanced, retarded and dual propagators, which is the basis of both, the duality 
theorem as a duality to the FTT, as well as the extension of the method to higher loop orders. The case of 
unstable particles with complex masses has been discussed in detail in Ref. [il2il . and we do not consider 
this possibility in the current paper. 

Besides the customary Feynman propagators Gp{qi), we also encounter advanced, GA{qi), and re- 
tarded, Gii{qi), propagators, defined by: 

qf -mj- iO qi^o ' ^ ' qf - mj + iO qt^ ' 

The Feynman, advanced, and retarded propagators only differ in the position of the particle poles in 
the complex plane. Using q]^ = qf^ — q^, we therefore find the poles of the Feynman and advanced 
propagators in the complex plane of the variable Qj o at: 



[GFiq^)^' = =^ %o = ±Vq? - mf - lO and [GAiq:)r' = =^ g.,o - ±^Jcd - mf+zO . (6) 

Thus, the pole with positive/negative energy of the Feynman propagator is slightly displaced below/above 
the real axis, while both poles of the advanced/retarded propagator, independently of the sign of the 
energy, are slightly displaced above/below the real axis (cf. Fig.[2l). We further define 

6{q,) ^ 2me{q,,o) S{q^ - m^) = 2m6+{q^ - m^) , (7) 

where the subscript + of 6+ refers to the on-shell mode with positive definite energy, gj o > 0. Hence, 
the phase-space integral of a physical particle with momentum g^, i.e., an on-shell particle with positive- 
definite energy, qf = mf, qi^ > 0, reads: 



d'^q 



^,0{%o)S{qt-mi) ■■■^ I 5{q,) ■■■ . (8) 



We continue by shortly recalling the duality theorem at one-loop order, which was derived in Ref. lfT2l . 
For a detailed discussion of all definitions and steps, as well as subtleties related to them, we refer the 
reader to this paper. In order to derive the duality theorem, one directly applies the residue theorem to the 
computation of L(^^(Pi,P2, • • • ,Pn) inEq. ©: Each of the Feynman propagators Gpiqi) has single poles 
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Figure 2: Location of the particle poles of the Feynman (left) and advanced (right) propagators Gpiqi 
and GA^qi) in the complex plane of the variable qi^. 



in both the upper and lower half-planes of the complex variable ii^. Since the integrand is convergent 
when £10 -^ cxo, by closing the contour at 00 in the lower half-plane and applying the Cauchy theorem, 
the one-loop integral becomes the sum of N contributions, each of them obtained by evaluating the loop 
integral at the residues of the poles with negative imaginary part belonging to the propagators Gpiqi)- 
The calculation of the residue of Gp{qi) gives 



Res[GF{qi)]im{q,,o)<o = / dii^oS+iqi - raf) 



(9) 



with 5+(gf — m^) defined in Eq. (|7]). This result shows that considering the residue of the Feynman 
propagator of the internal line with momentum q-i is equivalent to cutting that line by including the 
corresponding on-shell propagator (5+(gf — mf). The propagators Gp^qf), with j 7^ i, are not singular 
at the value of the pole of Gpiqi) and can therefore be directly evaluated at this point, yielding to 



where 



n ^^(9.0 



Goiquqi) 



Im(g,,o)<0 0^^ 



is the so-called dual propagator, as defined in Ref. [1121 . with 77 a future-like vector. 



(10) 



(11) 



(12) 



i.e., a d-dimensional vector that can be either light-like [rf = 0) or time-like (77^ > 0) with positive 
definite energy t/q. 

Collecting the results from Eq. ^ and Eq. (flOl) . the tree-loop duality theorem at one-loop, [[T2l . 
takes the final form 



L^^\puP2,---,Pn) 



E / ^{^^)\{GD{q^■.q,). 

"'^l 7 = 1 



(13) 






Contrary to the FTT, lfT4l [T5ll , Eq. (fT3l) contains only single-cut integrals. Multiple-cut integrals, like 
those that appear in the FTT, are absent thanks to modifying the original +iO prescription of the uncut 



Feynman propagators in Eq. (fTOl ) by the new prescription — iO r]{qj — Qi), which is named the 'dual' iO 
prescription or, briefly, the r] prescription. This is the main result of Ref. [12J. The dual iO prescription 
arises from the fact that the original Feynman propagator Gpi^j) is evaluated at the complex value of the 
loop momentum £i, which is determined by the location of the pole at (/^^ — m^ + iO = 0. The iO depen- 
dence of the pole of Gpiqi) modifies the iO dependence in the Feynman propagator Gpiqj), leading to 
the total dependence as given by the dual iO prescription. The presence of the vector 77^ is a consequence 
of using the residue theorem and the fact that the residues at each of the poles are not Lorentz-invariant 
quantities. The Lorentz-invariance of the loop integral is recovered after summing over all the residues. 
Furthermore, in the one-loop case, the momentum difference r]{qj —Qi) is independent of the integration 
momentum ii, and only depends on the momenta of the external legs (cf. Eq. (O). 

We now rederive the one-loop duality theorem Eq. (fT3l) by exploiting the relationship between Feyn- 
man, advanced and dual propagators. This will prove to be useful when extending the duality theorem 
to higher orders. Using the elementary identity 

PV(- ) TirrS{x) , (14) 



x±iO \x 

where P V denotes the principal-value prescription, we can transform one kind of propagators into the 
other: 

GAiQi) = Gpiqi) + 5 {qi) , Gniqi) = Gpiqi) + 5 {-qi) , ^^(-gi) = G^iqi) . (15) 

Dual and Feynman propagators are related through, [12], 

^(?i) Go^qu Qj) = ^(g*) Gpiqj) + {qj - q,) '5 {qj] 

with 9 (q) = 9{rjq). 



(16) 



In the following, we extend the definition of propagators of single momenta to combinations of 
propagators of sets of internal momenta: Let ak be any set of internal momenta with qi, qj G ak- We 
then define Feynman, advanced, retarded and dual propagator functions of this set ak in the following 
way: 

GF{A,R){ak) = Yl GF(A,R){qi) , Goiak) = ^ 6 (qi) JJ Goiquqj) ■ (17) 

By definition, G^icik) = ^(q'j)j when a^ = {{] and thus consists of a single four momentum. At 
one-loop order, ak is naturally given by all internal momenta of the diagram which depend on the single 
integration loop momentum ^1, afc = {1,2,..., N}. However, let us stress that a^ can in principle be 
any set of internal momenta. At higher order loops, e.g., several integration loop momenta are needed, 
and we can define several loop lines ak to label all the internal momenta (cf. Eq. (|T7l) ) where Eq. (fTTI) 
will be used for these loop lines or unifications of these. To simplify the notation, we also introduce 

Goi-oik) = ^ ^i-Qi) Yl ^oi-qu -qj) , GA{-ak) = Yl GA{-qi) = Gn^ak) , (18) 

where the sign in front of ak indicates that we have reversed the momentum flow of all the internal lines 
in ak- For Feynman propagators, moreover, Gp{—ak) = GF{ak). 



In analogy to Eq. (1151 ). the following relation holds for any set of internal momenta ak'- 

GAiak) = Gpiak) + Goiak) . (19) 

This is a non-trivial relation, considering Eq. (fT6l) . Note that individual terms in Gnictk) depend on the 
dual vector 77, but the sum over all terms contributing to Goictk) is independent of it. We leave the de- 
tailed proof by induction of Eq. (fT9l) for the Appendix lAl Eq. (fT9l) is our main result for a straightforward 
derivation of the duality theorem. 

Another crucial relation for the following is given by a formula that allows to express the dual func- 
tion of a set of momenta in terms of chosen subsets. Consider the following set /3jv = aiU...Uajv, where 
(3n- is the unification of various subsets ai. Solving for the dual part, Eq. (fT9l ) then has the following form: 



Goidi U ^2 U ... U ttAf) = 6*^(01 U a2 U ... U a^) — Gpicti U a2 U ... U ajy) . (20) 

We continue by using the multiplicativity of G^l/^Af) and Gf{Pn), as defined in Eq. (fTTI) . to obtain: 



N N 



G'i:)(ai U ^2 U ... U ajv) = TT ^^(aj) — TT ^^^(a 

1=1 i=l 

N 



N N 

i=\ i=\ 



Y. n ^^(«^i) n Gp{a,,). (21) 



The sum runs over all partitions of /3n into exactly two blocks /3]y and I3j^' with elements Oj, 
i G {1,...,A^}, where, contrary to the usual definition, we include the case: /3}^ = /S^, I3\^' = ^. 
This relation will be extensively used in the subsequent calculations. For the case of A^ = 2, e.g., where 

^2 = CKi U a2, we have: 

Gd{oli U 02) = GD{oii) GD{oi2) + Gd{oi.i) Gf{oi.2) + Gf{oli) G'z)(a2) • (22) 



Since in general relation (121]) holds for any constellation of basic elements «» which are sets of internal 
momenta, one can look at these expressions in different ways, depending on the given sets and subsets 
considered. If we define, for example, the basic subsets a^ to be given by single momenta g^, and since 
in that case Goiqi) = S (qi), Eq. (|2TI) then denotes a sum over all possible differing m-tuple cuts for 
the momenta in the set (3^, while the uncut propagators are Feynman propagators. These cuts start 
from single cuts up to the maximal number of cuts given by the term where all the propagators of the 
considered set are cut. 

Let us now return to the one-loop integral: Since advanced propagators have poles with positive 
imaginary part only, we do not enclose any singularity by closing the integration contour at 00 in the 
lower half-plane, and therefore a loop integral over advanced propagators vanishes. Thus, by using 
Eq. ^, we find 

= / Ga(«i) = / [Gpia,) + Goia,)] , (23) 

Jii Jii 

where ai as in Eq. ^ labels all internal momenta gj. The first term on the right-hand side of Eq. (|23l) . 
containing only Feynman propagators, is the original one-loop integral. Therefore, 

L«(pi,p2,...,P7v) = - / Gniai). (24) 




Figure 3: Momentum configuration of the two-loop N -point scalar integral. 

In this way, we directly obtain the duality relation between one-loop integrals and single-cut phase- 
space integrals and hence Eq. (I24l) can be interpreted as the application of the duality theorem to the 
given set of momenta ai. It obviously agrees with Eq. (fT3] ). Furthermore, by using Eq. (|2TI) in its 
refined form where the subsets Oj are given by the single momenta qi of the inner lines of the one-loop 
integral, we rederive the FTT at one-loop, namely the one-loop integral written in terms of multiple-cut 
contributions and Feynman propagators: 



L^^\pi,P2,---,Pn) 



E / n ^('^-) n ^H?: 

(2) "^ -t-l (1) (2) 



(25) 



^^(l)ua(^'=ai 



The sum runs over all partitions of ai as defined in Eq. (|2TI) . excluding the possibility to have a term 
with only Feynman propagators. The m-cut integral of the FTT is given by the sum of the contributions 
from all partitions of ai, with a{ containing precisely m elemental. 

The extension of the duality theorem and of the FTT from scalar loop integrals to full scattering am- 
plitudes in the case of unitary, local field theories and in the occurrence of real masses is straightforward 
and has been discussed in detail in Ref. [fT2l . 



3 Duality relation at two loops 



We now turn to the general two-loop master diagram, as presented in Fig.|3] Again, all external momenta 
Pi are taken as outgoing, and we have Pi J =pj+pi+i + . . .+pj, with momentum conservation pi at = 0. 
The label i of the external momenta is defined modulo A^, i.e., pN+i = Pi- Note, however, that one or 
both of the external momenta attached to the four-leg vertices might be absent: Pi = and/or pat = 0. In 
the two-loop case, unlike at the one-loop order, the number of external momenta might differ from the 



*If the number of space-time dimensions is d, then m is limited to be m < d; the terms with larger values of m vanish, 
since the corresponding number of delta functions in the integrand is larger than the number of integration variables. 



number of internal momenta. The loop momenta are ii and £2, which flow anti-clockwise and clockwise 
respectively. The momenta of the internal lines are denoted by qi and are explicitly given by 

{ii+PiA ,ieai 

i2+Pi,i-i ,iea2 (26) 

ii +£2+Pi,i-i ,i ^ as , 

where ak, with k = 1, 2, 3, are defined as the set of lines, propagators respectively, related to the mo- 
menta Qi, for the following ranges of i: 

ai = {0,l,...,r}, a2 = {r + l,r + 2,...,l}, a^ = {I + 1,1 + 2, ...,N} . (27) 

In the following, we will use a^ for denoting a set of indices or the set of the corresponding internal 
momenta synonymously. Furthermore, we will refer to these lines often simply as the "loop lines". 



We shall now extend the duality theorem to the two-loop case, by applying Eq. (|23l) iteratively. We 
consider first, in the most general form, a set of several loop lines ai to a^v depending on the same 
integration momentum £j, and find 

Gp{ai U 0:2 U ... U On) = — GDicti U a2 U ... U a^) , (28) 



which states the application of the duality theorem, Eq. (1231) . to the set of loop lines belonging to the 
same loop. Eq. (l28l) is the generalization of the duality theorem found at one-loop to a single loop 
of a multi-loop diagram. Each subsequent application of the duality theorem to another loop of the 
same diagram will introduce an extra single cut, and by applying the duality theorem as many times as 
the number of loops, a given multi-loop diagram will be opened to a tree-level diagram. The duality 
theorem, Eq. (l28l) . however, applies only to Feynman propagators, and a subset of the loop lines whose 
propagators are transformed into dual propagators by the application of the duality theorem to the first 
loop might also be part of the next loop (cf., e.g., the "middle" line belonging to 0:3 in Fig. [3]). The dual 
function of the unification of several subsets can be expressed in terms of dual and Feynman functions 
of the individual subsets by using Eq. (|2T|) (or Eq. ((22)) ') though, and we will use these expressions to 
transform part of the dual propagators into Feynman propagators, in order to apply the duality theorem 
to the second loop. 

We are now ready for extending the duality relation to two loops. Our starting point is the expression 
for the two-loop A^-leg scalar integral 

L'^^\pi,P2,...,Pn)= I / Gi.(ai U as U «3) , (29) 

where the momenta of the internal lines are specified in Eq. (|26l) and Eq. (1271 ). As stated before, we 
will apply the duality theorem sequentially to the two different loops associated with the integration 
momenta ii and £2. Starting with the first loop related to ii, and hence to the loop lines ai and 0:3, and 
using Eq. (|28l) . we obtain: 

L^^\pi,P2,...,Pn) = - GD{aiUa3)GF{a2). (30) 

Jii J 1.2 

We then use Eq. (|22|) for Goioii U a^), leading to 

L^^\pi,P2, ■ ■ ■ ,Pn) = - / {Gniai) Gnias) + Goiai) Gpias) + Gpiai) Gnias)} Gf{ol2) ■ 

(31) 



The first term of the integrand on the right-hand side of Eq. (|3T| ) is the product of two dual functions, and 
therefore already contains double cuts. We do not modify this term further. The second and third terms 
of Eq. (|3TI) contain one Gnicti) and hence single cuts only. We thus apply the duality theorem again, i.e., 
we use Eq. (l28l) for £2 in order to generate one more cut. A subtlety arises at this point since due to our 
choice of momentum flow, ai and 02, appearing in the third term of Eq. (|3T1) . flow in the opposite sense. 
Hence, in order to apply the duality theorem to the second loop we have to reverse the momentum flow 
of one of these two loop lines. We choose to change the direction of ai, namely g^ — )• —qi for i E ai. 
This change of momentum flow is denoted by a sign in front of a^. Thus, applying Eq. (128]) to the last 
two terms of Eq. (|3TI) . leads to 

n{-GDiai) GF{a2) Go^Ciz) + Goicii) G'z)(a2 U 03) + 6*0(03) Gd^-ch U 02)} . 

(32) 

This is the dual representation of the two-loop scalar integral as a function of double-cut integrals only, 
since all the terms of the integrand in Eq. (|32|) contain exactly two dual functions as defined in Eq. (fTT]) . 
The integrand in Eq. (|32|) can then be reinterpreted as the sum over tree-level diagrams integrated over a 
two-body phase-space. 

The integrand in Eq. (|32l) . however, contains several dual functions of two different loop lines, and 
hence dual propagators whose dual zO prescription might still depend on the integration momenta. This 
is the case for dual propagators Goifli] Qj) where each of the momenta q^ and qj belong to different 
loop lines. If both momenta belong to the same loop line the dependence on the integration momenta 
in 7]{qj — qi) obviously cancels, and the complex dual prescription is determined by external momenta 
only. The dual prescription ri{qj — q^) can thus, in some cases, change sign within the integration volume, 
therefore moving up or down the position of the poles in the complex plane. To avoid this, we should 
reexpress the dual representation of the two-loop scalar integral in Eq. (|32l ) in terms of dual functions of 
single loop lines. This transformation was unnecessary at one-loop because at the lowest order all the 
internal momenta depend on the same integration loop momenta; in other words, there is only a single 
loop line. 

Inserting Eq. (I22l) in Eq. (|32l ) and reordering some terms, we arrive at the following representation of 
the two-loop scalar integral 

L'^'^\pi,P2,---,Pn) 

n{GD{ai) G'D(tt2) Gp^d'i) + GD{-ai) G'F(a2) GdIos) + G*{ai) G'z)(a2) Gn^d'i)} , 
-. -2 

(33) 

where 

G*{au) = Gpiak) + Goiak) + Goi-ak) . (34) 

This is the second main result of this paper. In Eq. (|33] ). the iO prescription of all the dual propagators 
depends on external momenta only. Through Eq. (|34l) . however, Eq. (|33l) contains also triple cuts, given 
by the contributions with three Goictk)- The triple cuts are such that they split the two-loop diagram 
into two disconnected tree-level diagrams. By definition, however, the triple cuts are such that there is 
no more than one cut per loop line a^. Since there is only one loop line at one-loop, it is also clear why 
we did not generate disconnected graphs at this loop order. For a higher number of loops, we expect to 



find at least the same number of cuts as the number of loops, and topology dependent disconnected tree 
diagrams built by cutting up to all the loop lines a^. We explore this possibility at three loops in the next 
section. 

Note that using Eq. (fT9l ). G*{ak) can also be expressed as 

G*{ak) = GaM + GrM - Gpiak) , (35) 

which contains no cuts, although the imaginary prescription of the advanced and retarded propagators 
still depends on the integration loop momenta. 

Finally, let us remark that from Eq. (|33] ) we can obtain the FTT representation of the two-loop scalar 
integral. More precisely, we can write each Goifyk), for the lines k E {1,2, 3}, in terms of Feynman 
propagators and multiple cuts of their constituting momenta, by using Eq. (|2TI) with the basic subsets Oj 
given by single momenta. After regrouping some terms, we obtain 

l^'\pi,---,pn)= yi 1 1 [gfm n ^(9n) n ^(^^^) n ^^(^^3) 

fcG{l,2,3} 

+Gp(«2) n ^(-*i) n ^(^^^) n ^^(^^3) 
+Gi.(a3) n ^(^*i) n ^(^2) n '^^(^^3) 



+( n '^(*^)+ n '^(-'^n)) n '^(^-) n ^('^«) n ^^ 



iiGaj^ ii&a\ 12602 «3Ga3 i4£Oi Uofg Ucig 

(36) 



where the sums always run over the partitions of the index sets as defined and used in Eq. (|2TI) . We see 
that it consists of at least double cuts up to a multiple cut of all the internal momenta. 



4 Duality relation beyond two loops 

In this section, we will take a first look on the duality relation beyond the two-loop order by considering 
the master topologies at three loops, as represented in Fig. |4l We obtain explicit representations of the 
diagrams of Fig. |4]by using the iterative method described in the previous section. Although the dual 
representations obtained in this section are not unique, as the diagrams can be expressed in different ways 
in terms of dual and Feynman propagators depending on the choice of lines whose momentum flow is 
changed in the course of applying the duality theorem to each loop, we have followed a systematic way 
in order to minimize the number of terms. 

Diagrams Ht a) -He) are in a certain sense of the same type, and will be treated in the same way. We 
first cut these diagrams on the disjoint loops assigned to the lines {ai, 02 } and {0:3, 0:4}. Considering 
the basket ball diagram (a) of Fig. HI for example, this means: 

^balkct(Pi'P2, ...,Pn)= / / Goiai U 0:2) Goias U 04) . (37) 

Jii Jh Jh 
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(b) 



id) 



Figure 4: Master topologies of three-loop scalar integrals. Each internal line a^ can be dressed with an 
arbitrary number of external lines, which are not shown here. 

Remember that any dual function GD{oLk) of any set of momenta, and hence any application of the 
duality to a loop, contains at least one cut. Since we have a product of two expressions of the dual 
type, all terms in the expansion of this product in Eq. (|37] ) via Eq. (I22l) contain at least two cuts. 
Terms with triple and more cuts are already in their final form. They belong to the case where ei- 
ther all lines are dual, or one line is of the type "Feynman". However, the double-cut terms, as, e.g., 
(jD(ai) GF{ci2) Gpidz) Goicii), stemming from the combination of one Feynman propagator from 
each of the two different loops, form a third loop, which still consists only of Feynman propagators 
and hence still needs one more application of the duality theorem in order to generate the third cut: 

G D{ai) G F{ci2) G F^a^) G D{ai) ~^ - Goiai) G 0(021) as) G 0(0:4) ■ (38) 

UiJEiJh Jii Ji2 Je-z 

In the case of the zigzag diagram, Ub), or the ladder diagram, |3c), this third loop consists of one more 

internal loop line: line a^ in the former case and the lines a^ and a^ in the latter. Hence the third loop 

now consists of exactly one loop line from the first and second loop, and these additional loop lines. Due 

to the nature of the application of the duality through Eq. (|28l) and Eq. (|2TI) . we have to sum over all 

possibilities to build such sets of loop lines fulfilling these properties. Additionally, we have to assure 

that for each application of the duality the integration momentum runs in the same sense and hence 

change the momentum-flow direction for some chosen loop lines. Since each application of the duality 

generates one minus sign, we obtain the following general result for diagrams Ufa) -Stc): 




(3) 

^(a),(b),(c) 



(Pi,P2,...,Pn) 



mGD(ai U tta) GD(as U 04) Gf(P) 
-. -. -3 




Gd(oi2, 03, «4) ^'^(ai) + (^^(ai, tts, ^4) G'i?(a2) + Gd(oi.i, 0:2, 04) GF(ot-i) 

Gf(I^) 



+(?£>(«!, a2, 03) Gf(c^^ + Gfi(<^\-, 02, «3, «4) 



—Gd(o-\-, 03) Gd(o-2 U —ol\ U /?) — Gd(ol\, 0:4) Gd(o-2 U tta U /3) 



-GD(a2, as) GD(-ai U -0:4 U /3) - GD(a2, «4) (^^(-cn U as U /3) > , 



(39) 



where /3 = in the case of diagram Ht a) (with (7^(0) = 1), /3 = a^ for diagramH^b) and l3 = a^Ua^ in 
the case of the ladder diagramllfc). Note that for the ladder diagram, lines 0:5 and ag depend on the same 
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integration momentum and can therefore be considered as a single loop line as defined in this context. 
Hence, in this sense, the diagram naturally reduces to the zigzag-case, as long as the relative sense of 
momentum flow in these lines is correct and stays unchanged. For brevity, we defined the product of 
dual propagators as GD(ai, ...,aAr) = nj=i Goion), in contrast to Gd{oii U ... UaAr), given in Eq. (|2T]) . 
The dual representation of the three-loop scalar integral in Eq. (l39l) . contains mostly triple cuts. We have 
allowed for a single four-cut contribution in order to make this expression more symmetric, but this term 
can be rewritten in terms of triple-cut contributions with the help of Eq. (|22|) . 



If we expand all existing dual functions in Eq. (1391 ) in terms of dual functions of single loop lines by 
using Eq. (|2TI) . we obtain, e.g., for diagram Ufa): 

-^basket bi ' ^^2,- •• ,PAf) = - / / / < G'D(a2, Os, -"4) ^^(ai) + GDlttl , os, -"4) G'F(a2) 

Jh Je.2 Jh y 

+ G'z)( — «!, 02, ^4) Gpid^) + G'd( — «!, 02, «3) Gf^Ch) 

+ ^^(-ai, ^2, "3, «4) + ^^(ai, 02, as, -a^) + (^^(-a;!, q;2, "s, -"4) \ ■ (40) 



In this expression, the complex dual prescription of all the dual propagators depend on external momenta 
only, although at the price of generating disconnected tree diagrams. Similar results can be obtained for 
diagrams Hfb) andHfc). There are up to four cuts for diagramSfa) (cf. Eq. (1401)). five cuts for diagrams 
Sfb) and six cuts for diagram Sfc), although in this last case five cuts are enough if a^ U ag is considered 
as a single loop line. 

Also the Mercedes star diagrammed) can be expressed in terms of only three-cut contributions or in 
terms of three- up to six- cut contributions. However, due to the non-planar nature of this diagram, the 
way of obtaining its dual representation is slightly more involved, whereas the general idea as explained 
before stays the same. We achieve for the Mercedes star diagram the following dual representation: 




- (^^(ai, ^2, as) Gpioii, as, ae) + ^^(as U 0^4 U as) Gd{oli, 0:2) Gf{cx&) 
+ Gd(— «! U a4 U ae) Gd{(X2, as) ^^(as) + G'z)(— a2 U as U —Qq) Goicti, as) GF{a4) 
+ (^^(ai) [(?£)( as U a4) Gai^a^) G'ir(a2 U ag) — GD{a2 U as U a4 U ag) G£){ct5) 

— G£){a3 U a4) G'£)(— a2 U as U — ag)] 

+ G'£)(a2) [(?£)(— a 1 U ag) G£){a4) Gpi^a^, U as) — Goicti U as U as U — ag) Ga^ai) 

— Goi—cti U ag) G'z)(as U a4 U as)] 

+ GD(as)[G'z)(— a2 U as) GdI— ag) Gpiai U a^j — Goi—cti U — a2 U 04 U as) Gd(— ag) 
-G'D(-a2Ua5)G'i?(-ai Ua4Uag)] I . (41) 
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From the inspection of the three-loop case and from the general derivation of the method, it seems 
obvious how to extend it to even higher loop orders. Eq. (|2T1) can also be used to obtain the FTT 
representation of scalar integrals or scattering amplitudes at any loop order. 

The duality relation can be extended to evaluate not only scalar loop integrals, as discussed so far, but 
also complete Feynman diagrams. The extension of the one-loop duality relation from scalar integrals to 
Feynman diagrams was discussed in details in Ref. [ fT2| . This extension relies on the simple observation 
that the duality relation acts only on the Feynman propagators of the loop, leaving unchanged all the other 
factors in the Feynman diagram. This is valid in any unitary and local field theories. In spontaneously 
broken gauge theories, it holds in the 't Hooft-Feynman gauge and in the unitary gauge. In unbroken 
gauge theories, the duality relation is valid in the 't Hooft-Feynman gauge, and in physical gauges where 
the gauge vector n'^ is orthogonal to the dual vector tj^^, i.e., n • 77 = 0. This excludes gauges where n''' 
is time-like. At one-loop order, this choice of gauges avoids the appearance of extra unphysical gauge 
poles, which in other gauges (e.g. the time-like axial gauge) are also poles of the second order. Within 
the same choice of gauges, additional (unphysical) gauge poles are absent also at higher-loop level, and 
the duality relation can be straightforwardly extended from scalar integrals to Feynman diagrams. 

In Ref. [jl2] , it was also shown how the one-loop duality relation can directly be expressed at the 
level of full scattering amplitudes (or, more precisely, off-shell Green's functions). The derivation of 
the duality between one-loop and tree-level scattering amplitudes requires a detailed discussion of some 
issues related to tadpole and self-energy configurations. These (and related) issues become more delicate 
at higher-loop levels. We do not pursue further on this point in this paper, and we postpone detailed 
investigations to further studies. 



5 Conclusion and Outlook 

We have rederived the tree-loop duality theorem at one-loop order, which was introduced in Ref. [|T2l . 
in a way which is more suitable for extending it to higher loop orders. By iteratively applying the 
duality theorem, we have given explicit representations of the two- and three-loop scalar integrals. The 
method, however, is easily extendible to higher loop orders beyond three loops. In general, the dual 
representation of the loop integrals can be written as a sum of terms with exactly the same number of 
cuts as the number of loops, and in such a way that the loop diagram is opened to a tree-level diagram. 
However, this requires to deal with uncut propagators with complex dual iO prescription depending on 
the integration momenta, and thus with complex dual prescriptions that might change sign within the 
integration volume. Dual representations of the loop integrals with complex dual prescription depending 
only on the external momenta can be obtained at the cost of introducing extra cuts, which break the loop 
integrals into disconnected diagrams. This is a new feature of the duality theorem beyond one-loop, 
which does not happen at the lowest order. The number of extra cuts to be taken into account depends 
on the topology of the loop diagram. The maximal number of cuts agrees with the number of loop lines, 
and the cuts are such that it does not appear more than a single cut for each internal loop line. These 
general facts are true for the application of the duality relation to diagrams with an arbitrary number of 
loops. The results presented in this paper can also be used to obtain the FTT representation of diagrams 
at higher orders. 

The dual representations obtained in this paper are valid as far as only single poles are present when 
the residue theorem is applied. At one-loop order, the propagators of the gauge bosons might generate 
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unphysical poles, or even higher order poles. Those non-single poles can be avoided by a convenient 
choice of the gauge or of the dual vector [12]. At two- or higher loop orders, however, higher order 
poles might appear when diagrams with selfenergy insertions (nested or disjoint) are considered. At two 
loops this happens when two of the loop lines are made of single propagators, and no external momenta 
are attached to any of the two four-leg vertices of Fig. [3l At higher order loops, there are many more 
possible topologies showing this feature. Extending the duality theorem to this kind of diagrams requires 
to evaluate the contribution of the higher order poles, which depends on the topology of the diagram and 
on the nature of the internal propagators and on the form of the interaction vertices. Explicity studies of 
these loop diagrams are left to future investigations. 

Note added: After completion of this paper, a work ["201 dealing with similar topics appeared. The 
author uses retarded boundary conditions, and obtains some combinatorial factors weighting the different 
terms contributing to the loop-tree duality relation. We haved checked explicitly, in the two-loop case, 
that such combinatorial factors are the result of averaging over the different dual countertparts of the 
same loop integral obtained by permuting the loop lines Oj. Although many of the terms obtained in 
this way are either equivalent or can be related to each other after shifting the loop momenta, a larger 
amount of terms than by using a single dual counterpart are needed to be summed up for the same loop 
integral or Feynman diagram. By using Eq. (l55l) or Eq. (l57l) dual propagators can be expressed in terms 
of advanced or retarded propagators in a straightforward way, leading to an equivalent loop-tree duality 
relation as presented in Ref. [20] . In our paper we have not followed this procedure, since a main feature 
of our duality relation is that the 'zO' prescription of the dual propagators depends only on the external 
momenta. 
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A Derivation of some algebraic relations 

In this Appendix, we prove by induction several algebraic relations that have been used in the text. The 
basic ingredient of the proof is the following relation: 

^(Ai) ^(Ai + Aa) ... ^(Ai + A2 + ■ ■ ■ + A„_i) + cyclic perms. = 1 , (42) 

that holds for any set of n real variables Aj, with z = 1 , 2, . . . , tt,, such that 

n 

^A, = 0. (43) 

1=1 

Relation (|42l) was proven in Appendix B of Ref. [[T2]| . It applies, in particular, to Aj = 7]Pi, and follows 
from momentum conservation Yll=i Pi — 0' where pi are external momenta. In the following, we will 
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use Eq. (I42l) by setting Aj = r](qi — Qi+i) fori E {1, ..., n}, with (n + i) = i mod n. The four-momenta 
Qi are any arbitrary set of internal momenta, and the real variables Aj = ri(qi — q-i+i) might still depend on 
the loop momenta ii and £2- By construction, however, Eq. (|43l) is automatically fulfilled. Thus, Eq. (|42l) 
can also be written as: 

n n n n 

En^~(«^ - 1^^ = Y.U^(i^ - 9.) = 1 • (44) 

i=l j=l i=l J=l 

We start by deriving the following algebraic identity: 

Goioik) = GAiotk) - Gpiotk) , (45) 

where 

Goiak) = "^ HQi) n ^DiqhQj) , GA(F){o^k) = Yl GA{F){qi) ■ (46) 

as defined in Eq. (fTTj ). Remember that G/)(«fe) = ^ (?«) when a^ consists of a single four momentum q^. 
Using the identity GA{qi) = Gpiqi) + 5 {qi), the right-hand side of Eq. (|45l) . can also be written as 

g'aK)-g^k)= 5^ n ^(^*^) n ^^^(^^2), (47) 

where the sum runs over all possible partitions of a^ into exactly two subsets a\. ' and aj, ' . Additionally, 
we allow for the subset a^ to be empty, but the subset a\, always contains at least one element, this 



means that a term with only Feynman propagators is excluded. The first non trivial case for Eq. (|45l) 
occurs for a^ = {1, 2}. By using Eq. (fT6l) . and 9 (gi — ^2) + ^ {Q2 — Qi) = 1^ it is straightforward to 
prove that 

?(gi) Goiqi] q2) + ^(^2) GD{q2; Qi) = ?(gi) ^^.(^2) + ^(92) Gpiqi) + 6 {qi) 6 {q2) . (48) 

Let us now assume that relation (1451) is correct for A^ loop momenta, a^ = {1, . . . , A^}, and show 
that it is valid for a^^^ = a^ U {A^ + 1}. We replace the dual propagators Gc,{qi; qN+i) appearing in 
GDioi^~^^) by using again Eq. (fT6l) . We obtain: 



N+l N+1 



GdK^') = 5^(^(g.) n ^^(^*'^j) = ^^(«'^)^^(^^ 



i=l j=l 



N N N 



(iV iV iV \ 

Y, {qN+1 - qi) ^{qi) JJ ^oiqu qj) + H GniqN+i, q^) ) • (49) 



J=l J = 



Assuming that Eq. (|45|) and Eq. (|47l) are valid for A^ elements, the first term in the right-hand side of 
Eq. (|49l ). which is proportional to GpiqN+i), becomes 

Goioik) GpiqN+i) = [GA{ak) - Gpiak)] Gp^qN+i) ■ (50) 
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For the remaining terms in Eq. (|49l ). which are proportional to 6 (qN+i), we apply again Eq. (fT6l ) to all 
dual propagators. After some algebra, we find: 



GD{ar)-GD{ak)GF{qN+^) = 

5(gjv+i)(G^(«fc)+ Yl 0(«l'^U{iV + l}) H 6{q,,) H Gpiq,,)], (51) 

with 

e(a«U{iV + l})= 5^ f n HQn-Q^.)y (52) 

«27^n 

where the sum runs over all possible products of 9 functions that can be constructed with the four- 
momenta in the set a^ ' and qN+i- Obviously, from Eq. (|44l) . 

G(4'^U{iV +!}) = !, (53) 

for all possible partitions of ak into a^ ' and al '. Collecting the results from Eq. (l50l) . Eq. (ISTI) . and 
Eq. (l53l) . we finally obtain 

Goiak^^^) = GD{oik)GF{qN+i) + 5{qN+i){GF{ak) + GD{ak)) 

= GA{ak) GA{qN+i) - Gpiak) Gp^qN+i) , (54) 

as we wanted to demonstrate. 

Another useful relation of advanced and Feynman propagators is the following: 

gdK) = e ( n ^^(^^•) ) ^(^^) ( n ^^(^') ) • (55) 



We do not attempt to present a detailed proof here, as Eq. (1551) can be straightforwardly derived by 
reordering some terms from Eq. (l47l) . Analogously, we also have 



Goi-oik) = GR^ak) - Gpiak) , (56) 

and 

Goi-ak) = E ( n ^^('?^) ) ^(-^^) ( n ^^(*) ) ' (57) 

which also hold straightforward by changing the momentum flow of the four-momenta from ak in 
Eq. (l45l) and Eq. (l55l) and taking into account that Gpi—qi) = Gpiqi) and GA{—qi) = Gji{qi). 



B Massless sunrise two-loop two-point function 

We consider, as explicit example of the application of the duality relation at two-loops, the massless 
sunrise two-loop two-point function (Fig. [5]). From Eq. (l33l) . the dual representation of the sunrise 
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Figure 5: Sunrise two-loop two-point function. 
two-loop scalar integral is given by 

L^'\PUP2) = j j {^(?i) ^fe) GF{q,)+~6{-q^) GF{q2)~H(l^) + G*{qi)~5{q2) Jfe)} • (58) 
After replacing G*{qi) = Gpiqi) + S (qi) + S (— gi), and shifting some momenta, we obtain 

L^'\PI,P2) = I [ S{ii)S{i2) {GF{h + i2+Pl)+GF{il+i2-pi)+GF{ii-i2-pi) 

+S{h + i2+Pi) + S{ii + i2-Pi) ] . (59) 



For the integration of the first loop momentum £i, we use the basic integrals already calculated in 
Ref. El: 



Hi 

and 



/ 5 (i,) Gpiii + k) = dr [k^ + ^O] ' [l + Oik^) O^-k^) [e"^^' - l)] , (60) 



f 6 (£i) 6 {ii + k)=dr [P + iO] " 9{-P) (e^^"^ - l) , (61) 

where 

, _ cr 1 1 _ r(i + 6)r^(i-6) 

^ 2 e(l - 2e) cos(7re) ' ""^ (47r)2-^r(l - 2e) ' ^ ^ 

which we have reexpressed in a more suitable way in terms of [k"^ + iO] ^ Applying Eq. (l60l) and 
Eq. (ED to Eq. §^, we find 

L^'\pi,P2) = dr J^ Sih) I [{i2+Pi? + tOr (e^'"^ + l) 

+ [{i2-Pif + tO]-' [e^'"^ - 0{{i2-Pi?)em-Pi)o) (e^'^^ - 1)] ) • (63) 



The new phase-space integrals that we have to evaluate now for the integration over the second 
loop momentum £2 are then quite similar to those already encountered at one-loop. The calculation is 
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elementary, and we obtain 

dr [ 5ii2) [(^2 + A:)' + ^0] 
J19. 



'1-2 



sm(37rej '- \ /j 



and 



(64) 



dr / 5(4) [ii2 + ky + ^0^'eiii, + ky)eiii2 + k)o) = 

G2 ^^^ i-k' - i^r^^ [e{-k') - e{e) 9iko) e-^'-^] , (65) 

r - r(-i + 2e)r(i-6)^ 

^' - (4.r)4-2^r(3-3e) ' ^^^^ 



where 



Applying Eq. (1641) and Eq. (1651) to Eq. (1631) . and summing up all the theta functions, we finally get 

L^'\pl,P2) = -G2{-pl-^0Y-'', (67) 

which is the well-known result for the massless sunrise two-loop two-point function. 
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